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Abstract 

Asteroids and comets dissipate energy when they rotate about any axis 
different from the axis of the maximal moment of inertia. We show that 
the most efficient internal relaxation happens at twice the frequency of the 
body's precession. Therefore earlier estimates that ignore the double fre- 
quency input underestimate the internal relaxation in asteroids and comets. 
We suggest that the Earth's seismological data may poorly represent the 
acoustic properties of asteroids and comets as internal relaxation increases 
in the presence of moisture. At the same time, owing to the non-linearity 
of inelastic relaxation, small angle nutations can persist for very long time 
spans, but our ability to detect such precessions is limited by the resolution 
of the radar- generated images. Wobbling may provide valuable information 
on the composition and structure of asteroids and on their recent history 
of external impacts. 
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I. INTRODUCTION 



Comet P/Halley exhibits a very complex rotational motion (e.g., Peale and Lissauer 
1989), which is attributed to its rotation about an axis that does not coincide with the 
axis of the major inertia. Some asteroids also wobble like, for example, 4179 Toutatis 
(Ostro et al. 1993, Harris 1994, Ostro et al. 1995, Hudson and Ostro 1995, Scheeres et 
al. 1998, Ostro et al. 1999). 

A precessing body, as opposed to one steadily rotating about its axis of maximal 
inertia, is subjected to alternating stresses. These stresses deform the body, and cause 
energy dissipation due to inelastic effects. Naturally, this cannot change the angular 
momentum of a freely rotating body. Therefore the body tends to relax to the state of 
minimum energy, namely to the state of rotation with the maximum-inertia axis parallel 
to the angular momentum vector. External impacts may, however, excite wobble: if the 
body is observed precessing, this means that it has been subjected to impacts within the 
characteristic time of inelastic relaxation. Tidal interaction may be an additional source 
of wobble excitation. This source may be of a special importance for planet-crossing 
asteroids and comets (Black et al 1999). Jetting may be another reason for cometary 
precession. Generally, spinning bodies may start wobbling when they change their inertial 
axes through the loss of material (craters or outgassing). It is the balance between the 
impacts and tidal interactions, on the one hand, and the inelastic relaxation, on the other 
hand, that determines the dynamics of the body's rotation. 

Therefore a study of the rotation of asteroids and comets may provide valuable infor- 
mation about their recent history of external impacts (and tidal interactions), provided 
that we quantitatively understand the inelastic-dissipation process. 

An important study of inelastic relaxation in a wobbling asteroid has been performed 
by Burns & Safronov (1973). In their treatment they decomposed the complex pattern of 
body deformations into bending and bulge; calculated deformations and the elastic energy 
associated with the deformations, and calculated the energy-dissipation rate, using the 
material quality factor (so-called Q factor). 

Our results differ from those in Burns and Safronov (1973) because, first, the dissi- 
pation at the double-frequency was missed in the latter study (which in fact provides 
the leading inpunt in the dissipation process)^. Second, we shall try to be more exact 
in our quest for the probable values of the quality factor Q . Burns and Safronov in 
their article quoted the standard seismological data ( 10 2 < Q < 10 3 ), and eventually 
used Q ~ 300 as a more-or-less-acceptable average. Their listed range of values of 
Q is indeed typical for almost all terrestial rocks under the conditions natural for the 
lithosphere: at temperatures varying from the room temperature up to several hundred 
Celsius; at pressures about several MPa ; at frequencies from 10 _1 Hz through several 
kHz ; and most important, for rocks with traces of moisture. None of these conditions 
are believed to hold for tumbling asteroids and comets, and therefore we shall not be able 



Another attempt of calculating the stresses emerging in a precessing body, presented in a 
long-standing article (Purcell 1979), also omitted the double- frequency contribution to the stress 
tensor. 
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to borrow data from the geophysical literature. 

Uncertainties in the approach presented in the previous literature motivate our search 
for more rigorous methods, as well as for more appropriate values for the quality factor 
in order to provide more accurate estimates for the inelastic-relaxation time. 

This is the second paper in the series. In the first (Lazarian & Efroimsky 1999) 
we calculated the internal relaxation in a precessing dust grain. In what follows we shall 
remind the reader basic facts about solid body rotation (section 2). We shall calculate the 
stresses caused by the precession, and the rates of internal relaxation (all these technical 
calculations are in the Appendix). Then we shall compare our expression for the relaxation 
rate with that suggested by earlier researchers (section 3). Section 4 discusses whether 
our results are applicable not only to oblate but also to prolate bodies. A discussion of 
acoustical properties of materials is given in section 5. A particular example of wobble 
(asteroid 4179 Toutatis) is discussed in section 6, and the conclusions are summarised in 
section 7. 



We are interested in the free rotation of an oblate symmetric body. Its principal 
moments of inertia will be denoted by /«. One may assume, without loss of generality, 
that 



The body's inertial angular velocity will be denoted by f2, and its precession rate will 
be called uj. The body-frame-related coordinate system is naturally associated with the 
three principal axes of inertia: 1, 2, and 3, with coordinates denoted and unit 

vectors ei, e 2 , e 3 . The body- frame-related components of fl will be denoted as ^1,2,3 • 

The other, inertial, frame (X, Y, Z), with unit vectors ex, ey, e z , is chosen so that its 
Z axis is parallel to the (conserved) angular momentum J, and its origin coincides with 
that of the body frame (i.e., is at the body's centre of mass). Coordinates with respect 
to the inertial frame are denoted by the same capital letters as its axes: X, Y, and Z. 

We shall be interested on 9, the rate of the maximum-inertia axis' approach to the 
direction of angular momentum J. To achieve this goal, one has to know the rate of 
energy losses caused by the inelastic deformation. To calculate the deformation, we have 
to know the acceleration experienced by a particle located inside the body at a point (x, 
y, z). Note that we address the inertial acceleration, i.e., that with respect to the inertial 
frame (X, Y, Z), but we express it in terms of coordinates x, y and z of the body frame 



The fast processes (revolution and precession of a symmetric oblate body) are de- 
scribed by the Euler equations whose solution, ignoring any slow relaxation, reads (Fowles 
and Cassiday 1986, Section 9.5; Landau and Lifshitz 1976): 



II. SOLID-BODY ROTATION 



h > h 




(1,2,3). 



Oi = f2j_cosa;t 



f2 2 = sin out , f2 3 = const 



(2.2) 



where 
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fi_i_ = fi sin a , fi 3 = fi cos a 



(2.3) 



and 



fi±/fi 3 



tan a 



h tan # 



(2.4) 



Here 9 and a are the angles made by the maximal- inertia axis with J and fi . The 
angular velocity ft nutates around the principal axis 3 at a constant angular velocity 



u = (h- l)fi 



3- 



h = h/I . 



(2.5) 



However from the point of view of the inertial observer it is rather axis 3 that wobbles 
about J. Thus the angle 9 between J and axis 3 is constant^, as long as the energy is 
conserved. 

The components of ft are connected with the absolute value of the angular momentum: 



fi 3 = — - = — cos 9 



fi 



— h sin 9 

h 



(2.6) 



As ft preceses with rate u , its components bear a time-dependence in the forms 
of sin cut and cos cut . Since the centripetal component of the acceleration is quadratic 
in ft , a double frequency must unavoidably emerge in the expression for acceleration. 
Unfortunately, this circumstance has been missed in the literature hitherto, and all the 
authors have been considering dissipation at the principal frequency uj solely. In our 
recent article (Lazarian and Efroimsky 1999) we presented a comprehensive calculation 
of the acceleration. The time-dependent and time-independent components of the accel- 
eration give birth to time-dependent and time-independent components of the stresses, 
correspondingly. The stresses and strains arising at the double frequency considerably 
increase the elastic energy associated with the vibrations. This leads to a higher rate of 
energy dissipation in the body, and therefore to a much higher rate of relaxation. Calcu- 
lation of the alignment rate, for an oblate body modelled by a prism of sizes 2a x 2a x 2c, 
(c < a), is presented in the Appendix. Here follows the final result: 



de/dt 



3 • 3n 

- sin 3 



63 (c/a) 4 cot 2 9 + 20 



[1 + (c/a) 



214 



a 2 fig p 
li Q 



where fi is the shear elastic modulus of the material, Q is its quality factor, and 

fin = i 



(2.7) 



(2.8) 



is a typical angular velocity. The above formula (|2.7|) shows that the major-inertia axis 
slows down its alignment for vanishing 9 , which looks reasonable. Formula (|2.7j) differs 



2 The rate uj of precession is of the order of except in the case of h — > 1 , or in a very 
special case of ft and J being orthogonal or almost orthogonal to the maximal-inertia axis 3. 
Hence one may call the rotation and precession "fast motions" , implying that the relaxation is 
slow: 6 uj. It is in this sense that we assume 9 is constant 
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by a factor of 2 from the appropriate formula in (Lazarian and Efroimsky 1998) because 
of an error in our preceding article (see Appendix for details). 

It would be natural also to expect that the rate of alignment decays to zero for c/a 
approaching unity. (For c = a , the body simply lacks a maximum- inertia axis.) One 
may as well expect to see not only a "slow finish" but also a "slow start" (so that dd/dt 
vanish for 9 — > 7r/2 ): the maximum-inertia axis must be hesitant as to whether to start 



aligning along or opposite to the angular momentum. Still, if we look at ( |2.7|) , it will 
appear that dd/dt remains nonvanishing for c/a approaching unity, and that the major 
axis leaves the position 9 = n/2 at a finite rate. Recall, however, that all the above 
machinery works only insofar as the rotation and precession are fast, while the alignment 
is slow: 0<w . 

The time needed for the maximal-inertia axis to shift toward alignment with the 
angular momentum is: 

[Oo d9 

1 = ~J S d9/dt ■ (2,9) 

where 9q is the initial angle (9q < tt/2), while a small but finite S is introduced to avoid the 
"slow-finish" divergency. Each pair of values of 9q and 5 will give birth to one or another 
numerical factor of order unity, in the expression for the relaxation time. Calculations 
presented in Appendix B show that t is not very sensitive to the choice of angle 9q as long 
as this angle is not too small. (This weak dependence upon the initial angle is natural 
since the divergence emerges at small angles in the end of the motion.) A particular choice 
of 5 must be based on one's capability to recognise the precession by observational means. 
In the case of ground-based photometric experiments, one studies the amplitude of the 
lightcurve variation (which in the first approximation is proportional to the variation in 
the cross-sectional area of the body). A typical observational accuracy is 0.01 mag; that 
is, deviations from one rotation to the next less than 0. 01 mag cannot be confirmed to be 
real. From these assumptions, we arrive at some rough estimate of the minimal half-angle 
of the precession cone: 5 ~ 10° . Radar can produce images with resolution as fine as a 
decameter, in an absolute reference frame, and can discern the dimensions and spin states 
of asteroids in detail. In principle, with an optimum data set, an observation may reveal 
wobble with precession-cone half-angles of several degrees^. For most radars "several" 
means: about five degrees^. To be on the safe side, we shall take 

5 = 6° , (2.10) 

though, with the NEAR spacecraft scheduled to rendezvous with (433) Eros in the next 
year, nutations of only a degree or less may become detectable. 

By plugging ( |2.7| ) into ( |2.9|) one can get the relaxation time t as a function of 9q and 
c/a . A simple computation shows that t is not particularly sensitive to the oblateness 



3 Steven Ostro, private communication 1999 
4 Scott Hudson, private communication 1999 
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FIGURES 




FIG. 1. Comparison of the approximation (2.11) with the exact results for the relaxation 
time t. The exact calculation of t (scaled by a factor of pQ/a 2 Q,Qp ) was performed by 
combining formulae (2.9) and (2.7) for various realistic values of the parameter c/a . The final 
value of the angle 6 was assumed to be 5 = 6°, while its initial value was varying from 5 
through 7r/2 . 



parameter c/a when values of this parameter fall within the interval 0.5 - 0.9 (see Fig.l). 
As we can see from the plot, the dependence t(9o) may be approximated, for c/a = 0.5 
- 0.9, by the following expression that does not bear any dependence upon c/a : 

t ~ -2.8O0g + 7.84 # - 0.79 (2.11) 

This is a rough approximation. For the above value of 5 , it gives an error of about 
100 % for 9 ~ 5 and an error of about 10 % for 9 > ir/4. Nevertheless it is more 
than acceptable for our purposes, for it enables one to get an idea as to what a typical 
value of the relaxation time could be. 



III. COMPARISON WITH AN EARLIER STUDY 



We shall compare our result for the typical time of alignment with the estimate ob- 
tained in Burns & Safronov (1973): 

the numerical factor A being about 100 for bodies of small oblateness, and about 10 
for bodies of very irregular shapes. In astrophysics we normally encounter the case of 
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small oblateness (except perhaps some exotic species of cosmic-dust grains, called carbon 
flakes), so let us assume 0.5 ~ c/a ~ 0.9 . Our formula will give: 

ttZ result) « (1 - 2) f OT 9 « (2 - 3) 5 

*(oT™*t) « (3 - 4) /or £ « vr/4 (3.2) 



w (4 - 5) /or £ £ vr/2 



J (pur result) ~ I 41 p a 2 

while the aforequoted formula of Burns and Safronov will read: 

t* Q 
p a 2 fig 



We see that in the preceding study the relaxation time was overestimated, for small 
initial angles, by two orders. In other words, the effectiveness of the inelastic-relaxation 
process was underestimated by two orders. Even when the initial angle is not too small 
( 6 ~ 7r/4 ), the underestimation will be by a factor of 30 or so|. The three main 
reasons for this underestimate are the following. First, it is the contribution of the double- 
frequency mode missing in (Burns & Safronov 1973). In the sum (63(c/a) 4 cot 2 9 + 20) 
emerging in formula (|2.7|) , the term 63(c/a) 4 cot 2 9 appears due to the principal mode, 
while the term 20 appears due to the double-frequency mode. Thence, in the case of c 
close to a the double-frequency mode will give (after integration over 9) a considerable 
input, while for c <ti a (which is not a likely case for asteroids and comets, but may 
happen for cosmic dust) it provides an overwhelming contribution to the entire effect. The 
only case when the contribution of the second mode is irrelevant is the case of a small 
initial angle 9 . But in this case there exists another simple reason for Q3.1J ) to give a 
large error compared to our analysis: (|3.1| ) simply ignores any dependence upon the initial 
angle, and this is why it gives too long times for small angles. The third source of error 
in (ft.l|) is an assumption, accepted by Burns and Safronov, that the energy dissipation in 
the case of small oblateness is predominantly due to the bulge flexing. It seems that this 
simplifying assumption does not work well. 



IV. DYNAMICS OF PROLATE BODIES 



At first glance, the dynamics of a freely-spinning prolate body obeys the same prin- 
ciples as the dynamics of an oblate one: the axis of maximum inertia will tend to align 



5 In the case of body of a very irregular shape, the discrepancy between Burns and Safronov's 
result and our calculation will be less, because in this case Burns and Safronov suggest a smaller 
value for the numerical factor A in Q3.1| ). (See formula (23) in Burns & Safronov (1973), and a 
paragraph thereafter.) 
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itself parallel to the angular momentum. If we model a prolate body with a symmetric 
top, it will be once again convenient to choose it be a prism of dimensions 2a x 2a x 2c , 
though this time half-size c is larger than a , and thereby I3 = 1% > I\ . Then all our 
calculations formally remain in force, except that the right-hand side in ( |A19| ) changes its 
signF], provided we keep the notation 9 for the angle between J and the body-frame axis 
3 (parallel to dimension 2c ). This brings an illusion that axis 3 tends to stand orthogonal 
to J, which at first glance appears natural since axis 3 is now not the maximum- inertia 
but the minimum- inertia axis. 

Unfortunately, this extrapolation of the oblate-body description to a prolate-body case 
is of absolutely no practical relevance. The problem is that an arbitrarily small difference 
between ^3 and I2 will entail an entirely new character of precession (Synge & Griffith 
1959). For the first time this topic, in the context of asteroid and comet precession, was 
addressed by Black et al (1999). We are planning to dwell on the subject comprehensively 
in our next article. 

V. NUMBERS 

Values of the involved parameters may depend both upon the temperature of the 
body and the wobble frequency. Evidently, the temperature-, pressure- and frequency- 
caused variations of the density p are tiny and may be neglected. This way, we can 
use the (static) densities appropriate to the room temperature and pressure: p( slhcate ) « 
2500 kg m' 3 and p ( carbon ) » 2000 kg mT ? \ 

As for the adiabatic shear modulus fi , tables of physical quantities would provide 
its values at room temperature and atmospheric pressure, and for quasistatic regimes 
solely. As for the possible frequency-related effects in materials (the so-called ultrasonic 
attenuation), these become noticable only at frequencies higher than 10 8 Hz (see section 
17.7 in Nowick and Berry 1972). Another fortunate circumstance is that the pressure- 
dependence of the elastic moduli is known to be weak (Ahrens 1995). Besides, the elastic 
moduli of solids are known to be insensitive to temperature variations, as long as these 
variations are far enough from the melting point. The value of /x may increase by several 
percent when the temperature drops from room temperature to 10 K . Dislocations 
don't affect the elastic moduli either. Solute elements have very little effect on moduli 
in quantities up to a few percent. Besides, the moduli vary linearly with substitutional 
impurities (in which the atoms of the impurity replace those of the hosts). However 
hydrogen is not like that: it enters the interstices between the atoms of the host, and has 
marginal effect on the modulus. As for the role of the possible porosity, the elastic moduli 
scale as the square of the relative density. For porosities up to about 20 % , this is not 
of much relevance for our estimates^. 



6 For details see (Lazarian and Efroimsky 1999) 

7 We express deep thanks to Michael Aziz and Michael Ashby for their consultations on these 
topics. 
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According to Ryan and Blevins (1987), for both carbonaceous and silicate rocks one 
may take the shear-modulus value fx ~ 10 10 Pa . For the dirty ice Peale and Lissauer 
(1989) suggest, for Halley's comet, fx ~ 10 10 dyne/cm 2 = 10 9 Pa while Q < 100 . 

A proper choice of the values of the Q— factor for asteroids is to be the most prob- 
lematic subject. As well known from seismology, the Q— factor bears a pronounced 
dependence upon: the chemical composition, graining, frequency, temperature, and con- 
fining pressure. It is, above all, a steep function of the humidity which presumably affects 
the interaction between grains. The Q— factor is less sensitive to the porosity (unless the 
latter is very high); but it greatly depends upon the amount and structure of cracks, and 
generally upon the mechanical nature of the aggregate. Whether comets and asteroids 
are loose aggregates or solid chunks remains unknown. The currently available experi- 
mental data are controversial. On the one hand, it is a well-established fact that comets 
sometimes get shattered by tidal forces. Namely, Shoemaker-Levy 9 broke into 21 pieces 
on the perijove preceding the impact (Marsden 1993)0. It seems that at least some comets 
are loosely connected aggregates, though we are unsure if all comets are like this. 

As for the asteroids, the question is yet open, and the low density of asteroid 253 
Mathilde (about 1.2 g/cm 3 ) may be either interpreted in terms of the rubble-pile hy- 
pothesis (Harris 1998), or put down to Mathilde being perhaps mineralogically akin to 
low-density carbonaceous chondrites, or be explained by a very high porosity^. In our 
opinion, the sharply-defined craters on the surfaces of some asteroids witness against 
the application of the rubble-pile hypothesis to asteroids 0. And of course, we should 
mention here Vesta as a reliable example of an asteroid being a solid body of a structure 
common for terrestial planets: Hubble images of Vesta have revealed basaltic regions of 
solidified lava flows, as well as a deep impact basin exposing solidified mantle. Thus we 
would say that at least some asteroids are well connected solid chunks, though we are 



8 There is also a strong evidence in favour of Comet Brooks 2 having been rent into pieces 
by the jovian gravity in 1886 (Sekanin 1982). One may as well mention the disintegration of 
Comet West in 1976 (Melosh and Schenk 1993), though we would rather decline this argument 
because the comet was too close to the Sun, and was therefore warmed up. Catenae of aligned 
craters found on the lunar surface may become another possible evidence of comets being prone 
to shattering. The strongest argument in favour of the rubble-pile hypothesis comes from the 
recent discovery of such catenae on Callisto and Ganymede (Melosh and Schenk 1993). 

9 It should be emphasised at this point that in our opinion high porosity of a material does 
not necessarily imply this material being a rubble pile. There exist also an opposite viewpoint: 
according to Steven Ostro and Alan Harris (private communication 1999), high porosities (tens 
of percent) indicate that much of the body is unconsolidated. 

10 This is just our opinion. According to Steven Ostro's opinion (private communication 1999), 
the craters on Mathilde comparable to the object's radius can only be made in a rubble-pile 
asteroid. 
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uncertain whether this is true for all asteroids^ 

In the sequel, hence, we shall assume that the body is not a loosely connected aggregare 
but a solid rock (possibly porous but nevertheless solid and well-connected), and shall try 
to employ some knowledge available on attenuation in the terrestial and lunar crusts. 

We are in need of the values of the quality factors for silicate and carbonateous rocks. 
We need these at the temperatures from several K to dozens of K, , zero confining pressure, 
frequencies appropriate to asteroid precession (1CT 6 — 10~ 4 Hz), and (presumably!) 
complete lack of moisture. 

Much data on the behaviour of Q— factors is presented in the seismological literature. 
Almost all of these measurements have been made under high temperatures (from several 
hundred up to 1500 Celsius), high confining pressures (up to dozens of MPa), and un- 
avoidably in the presence of humidity. Moreover, the frequencies were typically within the 
range from dozens of kHz up to several MHz. Only a very limited number of measure- 
ments have been performed at room pressure and temperature, while no experiments at 
all have been made thus far with rocks at low temperatures (dozens of K ) . The informa- 
tion about the role of humidity is extremely limited. Worst of all, only a few experiments 
were made with rocks at the lowest seismological frequences (10~ 3 — 10 _1 Hz), and none 
at frequencies between 10~ 6 and 10~ 2 Hz, though some indirectly achieved data are 
available (see Burns 1977, Burns 1986, Lambeck 1980, and references therein). 

It was shown by Tittman et al. (1976) that the Q— factor of about 60 measured under 
ambient conditions on an as-received lunar basalt was progressively increased ultimately 
to about 3300 as a result of outgassing under hard vacuum. The latter number will be 
our starting point. The measurements were performed by Tittman et al. at 20 kHz 
frequency, room temperature and no confining pressure. How might we estimate the values 
of Q for the lunar basalt, appropriate to the lowest frequencies and temperatures? 

As for the frequency-dependence, it is a long-established fact (e.g., Jackson 1986, 
Karato 1998) that Q ~ uo a with a around 0.25. This dependence reliably holds for all 
rocks within a remarkably broad band of frequencies: from hundreds of kHz down to 
10~ 1 Hz. Very limited experimental data are available for frequencies down to l§~ z Hz, 
and none below this threshold. Keep in mind that a being close to 0.25 holds well only 
at temperatures of several hundred Celsius and higher, while at lower temperatures a 
typically decreases to 0.1 and less. 

As regards the temperature-dependence, there is no consensus on this point in the 
geological literature. Some authors (Jackson 1986) use a simple rule: 

Q ~ u a exp(A*/RT) , (5.1) 

A* being the apparent activation energy. A more refined treatment takes into account the 
interconnection between the frequency- and temperature-dependences^. Briefly speak- 
ing, since the quality factor is dimensionless, it must retain this property despite the 



Currently available radar and optical data establish that the 30-meter object 1998 KY26 
is monolithic, and suggest that at least a few larger objects are also monolithic (Ostro 1999). 

12 The authors are grateful to Shun-ichiro Karato who drew our attention to this connection. 
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exponential frequency- dependence. This may be achieved only in the case that Q is 
a function not of the frequency per se but of a dimensionless product of the frequency 
by the typical time of defect displacement. The latter exponentially depends upon the 
activation energy, so that the resulting dependence will read: 

Q ~ [u exp{A*/RT)] a , (5.2) 

where A* may vary from 150 - 200 kJ/mol (for dunite and polycristalline 
forsterite) up to 450 kJ/mol (for olivine). This interconnection between the frequency- 
and temperature-dependences tells us that whenever we lack a pronounced frequency- 
dependence, the temperature-dependence is absent too. It is known, for example (Bren- 
nan 1981) that at room temperature and pressure, at low frequencies (10 -3 — 1 Hz) 
the shear Q— factor is almost frequency-independent for granites and (except some spe- 
cific peak of attenuation, that makes Q increase twice) for basalts. It means that within 
this range of frequencies a is small (like 0.1, or so), and Q may be assumed almost 
temperature-independent too. 

Presumably, the shear Q-factor, reaching several thousand at 20 kHz , descends, in 
accordance with (|5.2p , to several hundred when the frequency decreases to several Hz. 
Within this band of frequencies, we should use the power a ~ 0.25 , as well known 
from seismology. When we go to lower frequencies (from several Hz to the desirable 
10~ 6 — 10~ 4 Hz), the Q-factor will descend at a slower pace: it will obey ( p.2| ) with 
a < 0.1 . Low values of a at low frequencies are mentioned in Lambeck (1980) and 
in Lambeck (1988)Q. The book by Lambeck containes much material on the Q-factor 
of the Earth. Unfortunately, we cannot employ the numbers that he suggests, because 
in his book the quality factor is defined for the Earth as a whole. Physically, there 
is a considerable difference between the Q-factors emerging in different circumstances, 
like for example, between the effective tidal Q-factorf^] and the Q-factor of the Chandler 
wobble). In regard to the latter, Lambeck (1988) refers, on page 552, to Okubo (1982) 
who suggested that for the Chandler wobble 50 < Q < 100 . Once again, this is a 
value for the Earth as a whole, with its viscous layers, etc. We cannot afford using these 
numbers for a fully solid asteroid. 

Brennan (1981) suggests for the shear Q— factor the following valuesff} Q^shear)^ ~ 

250 , Q^shear) ~ 500 . It would be tempting to borrow these values^, if not for one 
circumstance: as is well known, absorption of only several monolayers of a saturant may 



13 See also Knopoff (1963) where a very slow and smooth frequency-dependence of Q at low 
frequencies is pointed out. 

14 A comprehensive study of the effective tidal Q-factors of the planets was performed by Gol- 
dreich and Soter (1966) 

15 Brennan mentions the decrease of Q with humidity, but unfortunately does not explain how 
his specimens were dried. 

16 Recall that these data were obtained by Brennan for strain amplitudes within the linearity 
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dramatically decrease the quality factor. We have already mentioned this in respect to 
moisture, but the fact is that this holds also for some other saturantsQ. Since the asteroid 
material may be well saturated with hydrogen (and possibly with some other gases), its 
Q-factor may be much affected. 

It may be good to perform experiments, both on carbonaceous and silicaceous rocks, 
at low frequencies and temperatures, and with a variety of combinations of the possible 
saturants. These experiments should give us the values for both shear and bulk quality 
factors. The current lack of experimental data gives us no choice but to start with the 
value 3300 obtained by Tittman for thoroughly degassed basalts, and then to use formula 
( ^2|) . This will give us, at T = 20 K and lu = 10" 5 Hz : 

- 150 (5.3) 

This value of the shear Q-factor^ for granites and basalts differs from the one chosen in 
Burns and Safronov (1979) only by a factor of 2. For carbonaceous materials Q must be 
surely much less than that of silicates, due to weaker chemical bonds. So for carbonaceous 
rocks we shall choose the following upper boundary: 

Q(carb) < 10Q ^ ( 54 ) 

though this boundary is somewhat arbitrary, and is probably still too high. 

Terrestial geophysics has not yet given us a reliable handle on the acoustic proper- 
ties of asteroids and comets, though the necessary quality factors can be obtained by a 
modification of the existing testing technique. Namely, the frequencies must be within 
the interval 10 -6 — 10~ 4 Hz, the temperature must be from about K to dozens of 
K . The specimens must be well outgassed and the measurements must be performed 
under a high vacuum that would mimic the real interplanetary environment. It would be 
most important to study specimens exposed to a variety of saturants (hydrogen, first of 
all). The specimens must include both silicate and carbonaceous rocks, and artificially or 
naturally prepared samples of dirty ice. 

Twenty two years ago J. A. Burns (1977) wrote: "Further experimental and theoretical 
work for real materials is sorely needed if we are to ever trace the history of the natural 
satellites." This plea is even more relevant today, and the ramifications of such work will 
nowadays be even broader: these studies are to become our key not only to satellite but 
also to asteroid studies. 



range. Our case is exactly of this sort since the typical strain in a tumbling body will be about 
a/n ~ p$7 2 a 2 //i . For the size a ~ 1.5 x 10 4 m and frequency not exceeding 10 -4 Hz, this 
strain is less than 10~ 6 which is a critical threshold for linearity. 

17 like, for example, ethanol (Clark et al. 1980) 

18 The bulk Q-factor may differ from the shear one. In our case we have a sophisticated de- 
formation picture that includes both torsional and longitudinal displacements. Simply from 
looking at the expressions for stresses we can see that torsion will dominate. Hence the effective 
Q-factor must be close to the shear Q-factor value. 
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VI. PARTICULAR EXAMPLE: ASTEROID 4179 TOUTATIS 



The asteroid 4179 Toutatis is a slowly rotating body of size about 2 km, that was 
"imaged" in radar by Ostro et al. (1993), (1995), (1999), Hudson and Ostro 1995, Scheeres 
et al. 1998. It is of S-type, analogous to stony irons or ordinary chondrites but certainly 
not to carbonaceous chondrites (Ostro et al. 1999). It has a rotation period r = 
7.5 days (so that fi = 2vr/r « 9.7 x 10~ 6 Hz). 

Our analysis was aimed at oblate bodies and, rigorously speaking, it is unapplicable to 
prolate or triaxial rotators. Still, while a similar study for triaxial rotators in on the way, 
let us try to apply our formulae, as a zeroth approximation, to Toutatis. Suppose that 
some tidal interaction has forced this asteroid to precess with a small precession-cone half 
angle 6 . Let it be, for example, twice or thrice the minimally recognisable half-angle: 
#o = (2 — 3) 5 ~ 12° . Then Fig.l (or our formula ( |3.2p based on it) will give, for 
Q = 150 , /I = 10 11 dyne/cm 2 = 10 10 Pa , and p = 2000 kg/m 3 , c/a ~ 1/2 : 

t ~ 2.5 ^„ = 0.5 x 10 18 s w 1.6 x 10 10 years (6.1) 
pa 1 \l% 

which is 800 times less than the estimate that would come from Burns & Safronov's 
(1973) treatment. Factor 800 emerges, first, due to the difference by a factor of 40 
between our formula fl6.1| ) and Burns & Safronov's formula (|3.3D; and second, due to 



the difference in choice of the quality-factor and shear-modulus values. Our choice is: 
p,Q = 1.5 x 10 13 dyne/ cm 2 = 1.5 x 10 12 Pa , while Burns & Safronov suggested 
p,Q = 3 x 10 14 dyne j cm 2 = 3 x 10 13 Pa . We suggest values lower than in Burns & 
Safronov (1973) for the reasons explained above. 

Harris (1994) goes even further: he suggests to use data measured for Phobos: fiQ = 
5 x 10 12 dyne/cm 2 = 5 x 10 11 Pa which is almost two orders of magnitude less 
than that in Burns & Safronov (1979). Therefore if we adopt the values proposed by 
Harris (1994), our estimate ( |6.1|) will be further decreased by 3 times, and will equal to 
5 x 10 9 years , which is 10 3 times less than the result that would follow from Burns & 
Safronov's treatment. 

Even though Harris chose a lower value for p,Q than we did, our calculation will give 
a lower value for the relaxation time, due to the difference in numerical factors in our 
formula ( |3.2|) and Burns & Safronov's formula (|3.3| ). (Actually, this difference is even 
larger than it should be, because Harris, applying Burns Sz Safronov's formula, took the 
overall numerical factor to be not 10 but 30.) This would mean that the number of 
asteroids that are suspected of tumbling (see Fig. 1 in Harris (1994)) should be decreased. 

How certain is this conclusion? It follows from our discussion in the previous section 
that at least for solid asteroids the piQ factor obtained for Phobos in Yoder (1982) may be 
an underestimate. Observations of asteroid precession should provide us with insight into 
their mechanical properties. The comparison of these with the properties of laboratory- 
studied materials should provide us with more understanding of the composition and inner 
structure of asteroids. 

While studying wobbling, it is important to bear in mind two points. First of all, 
nutations at small angle can proceed for much longer times than we estimated above. This 
is the consequence of the "slow finish" condition that we discussed in section III. With 
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the improvement of the radar measurements, we expect to detect more asteroids wobbling 
at small angles. Therefore high-accuracy imaging of asteroids may reveal nutations above 
the upper curve in Fig. 1 in Harris (1994). Second, large-amplitude or/and rapid tumbling 
may be suppressed quicker than we estimated as the strains become nonlinear and the 
dissipation increases^. 

VII. CONCLUSIONS 

1. The dissipation at the double frequency of precession makes an importantm, some- 
times decisive, contribution into the process of inelastic relaxation. 

2. Small-amplitude wobbling can proceed for long periods of time, while the large- 
amplitude precession gets damped to smaller amplitudes much faster. This is a conse- 
quence of the dependence of the relaxation rate upon the angle of precession. In reality, 
the finite resolution of radar-generated images makes us take into consideration only wob- 
ble with precession-cone half angles no less than about 5 — 6° . (Though observations 
made by spacecrafts will, hopefully, improve the resolution up to one degree or so.) 

3. Neglection of the above two circumstances, along with acceptance of an assumption 
(unjustified, in our opinion) that the tidal-bulge flexing plays a decisive role in the inelastic 
relaxation, lead our predecessors to a great underestimate of the effectiveness of the 
inelastic dissipation mechanism. (By two orders, for small initial angles, and by a factor 
of 20 - 30 for large angles.) 

4. Solid asteroids are likely to have Q— factors about one hundred, or possibly even less, 
but further laboratory research is required. Earth seismological data may poorly represent 
acoustic properties of asteroids and comets as even a tiny trace of moisture (several 
monolayers) substantially alters the quality factor. Laboratory data on 10~ 6 — 10~ 4 Hz 
oscillations of de-moisturised silicate and carbonaceous rocks at temperatures within the 
range from several K to several dozens K are badly needed. 

5. Potentially, wobbling may provide valuable information on the composition and 
structure of asteroids, and on their recent history of external impacts. 

6. Since the inelastic dissipation turns out to be a far more effective process than 
believed previously, the number of asteroids expected to wobble with a large angle of the 
precession cone is expected to be lower (than the predictions made in (Harris 1994).) 
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exceeding the linearity threshold. 
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APPENDIX A: ALIGNMENT OF THE MAXIMAL-INERTIA AXIS TOWARDS 
THE ANGULAR-MOMENTUM VECTOR 



If we model the freely rotating oblate body by a prism of sizes 2a x 2a x 2c, (c < a) 
and density p , the time-dependent stresses will read (Lazarian and Efroimsky 1998): 



pQ] 



(x — a ) cos 2ut , a. 



P n 2 _ 



\y 2 - a 2 ) cos2ut , a zz = (Al) 



a 



P n2 



■HI 



Q 2 ± (x 2 + y - 2a 2 ) sin 2ut , 



(A2) 



/i (2 2 - c 2 ) + (2 - h) {x 2 - a 2 ) 



cos cut 



(A3) 



«2 



| fi ± fi 3 [ /i - c 2 ) + (2 - /i) (y 2 - a 2 ) 



sin cut . 



(A4) 



As expected, not only the principal frequency but also the second mode shows up in cr^. 
The moment of inertia I3 and the parameter h read: 



r 16 4 

^3 = y P a c , 



According to ( pT6| 



fi cos 6* 



h = 7 



1 + {c/af 



Qn h sin 8 



(A5) 



(A6) 



Q = J/I 3 being the typical angular velocity of a body. We shall compute the strain 
tensor in order to estimate the maximal elastic energy W stored in the body. Then, 
we shall estimate the dissipation rate as — uW/Q, Q being the quality factor of the 
material. 

The kinetic energy of an oblate spinning body reads, according to Q2.ip, (|2.2|), and 
^61): 



E, 



rot 



[i Ql + i 3 n 3 



sin 



cos 2 e 



j 



so that 



dE„ 



rot 



J 2 

— (h — 1) sin 9 cos 9 = to J sin 9 



d9 , 3 

The rotational energy changes via the inelastic dissipation: 

E rot = W , 

W standing for the elastic energy of the body. Then the rate of alignment is: 



(A7) 



(A8) 



(A9) 
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d6 ( dE rot \ dE rot . ! • /A in \ 

— (oj J sin v) W (A10) 



dt V dd J dt 
where 

W = WW + W'^) = -w^- - 2^^^ « ^ {WW + 2W/^} , (All) 

the quality factor being almost frequency-independent. In reality, it certainly does depend 
upon u, but the dependence is very smooth and may be neglected for frequencies differing 
by a factor of two. (It should be taken into account within frequency spans of several 
orders.) In (|A11|) Wtf , Wq w denote the amplitudes of the elastic energies associated with 



vibrations at the modes u and 2u, while W u = 2Wq and W 2uJ = 2Wq uj stand for 
the averages. In our preceding publication (Lazarian and Efroimsky 1998) we missed the 
coefficient 2 between Wq and W^"\ 

Now the question is: how to calculate + 2W^ 2u1 ^ ? At low temperatures the 

bodies manifest, for smalQ displacements (like, say, sound) no viscosity: urj ~ u( <C 
H ~ K . Hence the stress tensor will be approximated, to a high accuracy, by its elastic 
part: 

Tr a f 1 \ 1 

6ij = 5ij ~9~K + 1°^ ~~ 3 5ij Tra ) 2~Jl ' ^ A12 ^ 

K and // being the bulk and shear moduli. The elastic energy stored in a unit volume is 

dW/dV = icyrry = I{(!£-I) (Tra) 2 + cr^^j . (A13) 

Anticipating different rates of dissipation in the two modes, we split dW/ dV into parts: 

dW/dV = dW {u) /dV + dW^/dV (A14) 
According to ( |A11| ), what we need is the sum 

+ 2W {2ul) = J a dx J a dy J C dz [dW^/dV + 2 dW {2ul) /dv} (A15) 



To compute this sum, one has to plug the expressions for stresses ( |A1[) - ([AJ) into (|A13|) , 
then to split the result in accordance with (|A14| ), and to build up the sum under the 
integral in (|A15|) . These calculations (presented in (Lazarian and Efroimsky 1999)) yield: 



>2\ 2 



W M + 2W ™ = 3 2 2- 6 (63 (c/a) 4 tan^g + 20) _^c^_ f 1^ . ^ 



V [1 + (c/a) 



(A16) 



20 The term "small deformation" means not exceeding the elastic limit. According to Brennan 
(1981), such deformations cause strain amplitudes not exceeding (for most materials) (1 — 3) x 
1(T 6 
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where Q = J/J 3 is a typical angular velocity. Eventually, (|A10|) , ( |A11| ) and ( |A16|) will 
yield: 



dd/dt 



According to ( |A5| ) and (|2.8| 



3 2 2" 6 (63 (c/a) 4 tan" 2 # + 20) a" 2 <r ] 



[1 + (c/a) 2 ]' 



sin 3 (A17) 



i //s^r 2 



J 



- ft 3 p a 4 c 



(A18) 



Substitution of the latter in the former gives the final expression for the alignment rate: 

(A19) 



3 , n 63 (c/a 4 cot 2 9 + 20 a 2 ft 3 o 

= S in 3 ^ — ' , , xo1/1 

' 2 4 1 + (c/a) 24 pQ 
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